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Just as big bang nucleosyntliesis allows us to probe the expansion rate when the 
temperature of the universe was around 1 MeV, the measurement of gravity waves 
from electroweak scale first order phase transitions may allow us to probe the ex- 
pansion rate when the temperature of the universe was at the electroweak scale. We 
compute the simple transformation rule for the gravity wave spectrum under the 
scaling transformation of the Hubble expansion rate. We then apply this directly to 
the scenario of quintessence kination domination and show how gravity wave spectra 
would shift relative to LISA and BBO projected sensitivities. 



1. INTRODUCTION 



The detection of gravity waves (GWs) generated during a scalar sector's first order phase 
transition (PT) represents an interesting future possibility [Tl-[26|. During a first order PT, 
bubbles of true vacuum nucleate, stir up the cosmological fluid, and collide, producing GWs. 
First order PTs at the electroweak scale in particular have received much attention because 
of their possible connection to a well motivated electroweak baryogenesis scenario in which 
the baryon asymmetry is generated during an electroweak phase transition (EWPT) |27t[28]. 
Studies of this scenario are particularly timely given that Tevatron and the LHC are actively 
probing the Higgs sector responsible for electroweak symmetry breaking. If the Higgs boson 
or any degree of freedom that can be responsible for a first order PT at the electroweak scale 
is found in the ongoing experiments, future experiments may eventually be able to measure 
all of the parameters necessary to give an accurate prediction for the GWs. It is important to 
emphasize that even if the electroweak symmetry breaking is not a first order PT, a typical 
beyond the standard model scalar sector has multiple degrees of freedom, and some of these 
can undergo a first order PT at a temperature near the electroweak scale. 

Just as the relative isotope abundance measurements have led to a constraint on the 
expansion rate during big bang nucleosynthesis, the measurement of GWs may allow us 
to constrain any non-standard expansion rate during the time of EWPT. Several detailed 
computations of the gravity wave spectrum exist, and each give varying degrees of dependence 
on the Hubble expansion rate. However, to our knowledge, previous work does not sufficiently 
discuss the general dependence of the Hubble expansion rate to directly answer the following 
question: how would the observed gravity wave spectrum change if the Hubble expansion 
during the EWPT was changed from that of pure relativistic degrees of freedom? 

We compute a simple transformation rule for the gravitational wave spectrum in terms 



* Electronic address: 
^Electronic address: 



danielchung@wisc.edu 


pzhou@wisc.edu 





2 



of ^ = H^/Hiji where Hu is the expansion rate which assumes radiation domination and 
is the actual expansion rate: 



implies the following: 1) The peak frequency of the spectrum will shift from the standard 
scenario frequency fp as fp — )■ ^/p, and 2) the peak amplitude of the spectrum will be 
suppressed from the standard scenario amplitude Ap as Ap — )■ Ap/^"^. The intuition for the 
amplitude is that less source contributes to the gravitational wave at a typical spacetime point 
today because of the smaller intersection of the past null boundary with the approximately 
compact time support of the source. The intuition for the frequency shift is that all conformal 
symmetry breaking scales relevant for the observable frequency range is controlled by the 
Hubble expansion rate. 

We then apply this scaling relationship to the results of [T^ [T6| to compute how the grav- 
ity wave spectrum will transform due to the assumption of the existence of a quintessence 
kination dominated phase |29]435] . Such assumptions are interesting because as pointed out 
by [32j (related scenarios were also suggested before by [HP, the freeze-out abundance 
of thermal relics can be strongly enhanced in scenarios in which the energy density is domi- 
nated by the kinetic energy of the quintessence field (kination domination) during the time 
of freeze-out, but dilutes away by the time of big bang nucleosynthesis (BBN). Such kination 
dominated freeze-out scenarios are then consistent with standard cosmology and predict that 
the standard relic abundance computed from the parameters extracted from collider mea- 
surements will be mismatched from the relic abundance deduced by observational cosmology. 
Thus, this scenario has interesting implications for physics models with thermal dark matter 
candidates (e.g. models with low energy supersymmetry such as the minimal supersym- 
metric extension of the SM (MSSM), technicolor models, models with large/warped extra 
dimensions, or certain classes of little Higgs models), which will be probed at the LHC and 
other experiments in the foreseeable future (for collider implications of this class of models, 
see for example [37]). Furthermore, as pointed out by |38l l39]. large annihilation cross sec- 
tions compatible with the dark matter explanation of the excess positrons [101 can be 
compatible with the right thermal relic abundance since the effective boost factor coming 
from the kination scenario can easily be as large as 10^. Finally, since the measurement of 
CMB B-mode polarization can almost model independently falsify this scenario, this sce- 
nario can be nontrivially checked with a variety of probes. In particular, the EWPT gravity 
wave probe in conjunction with dark matter cosmology can represent a smoking gun probe 
of the scenario if the gravity wave is measurable [38j and colliders can eventually measure 
the requisite short-distance parameters with sufficient accuracy. 

The order of presentation is as follows. In the next section, we present the main analytic 
result of this paper. Sec. |3] focuses on checking explicit consistency of our result with the 
existing literature on explicit gravity wave spectrum computations. In Sec. |4[ we apply the 
transformation relation to the quintessential kination scenario and give plots showing how 
the transformed spectra look relative to the projected sensitivities of LISA and BBO. Sec|5] 
discusses all of the caveats associated with the analytic result. We then conclude with a 
summary. The appendices contain some of the details used throughout the paper. 

As far as conventions are concerned, we use the reduced Planck's constant Mp = 2.4 x 10^^ 
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GeV and also assume a flat FRW background metric ds^ — aF'{t){dt^ — 



2. GENERAL ANALYTIC ARGUMENTS 

In this section, we compute how the gravitational wave spectrum will transform under the 
situation that during the last part of the PT, the Hubble expansion rate is ^Hu where Hu 
is what the expansion rate would be in the "usual" radiation domination epoch. Because the 
main arguments rely only on the general form of the gravity wave equation and dimensional 
analysis, the transformation results will be very robust and nearly model independent. In 
a latter section, we check the consistency of the transformation rules with explicit model 
dependent computations in the literature. 

Consider the transverse-traceless perturbation hj^^ about the background metric: 

ds"^ = a\t) [dt"^ - {6ij + hj^)dx'dx^] . (2) 

Using the pseudotensor expression for the energy density in gravity waves, the energy density 
in GWs can be expressed as 

PGw-^{doh^;'dohj;'). (3) 

During a phase transition (PT) the energy density in the gravity wave can be written as 

PGw{x,t) ~ (^)' (|^[A,^.,„| d'x'Gretix;x')TUx')]^')\pT (4) 

where a* is the scale factor at the beginning of the PT, Gj-Q\^{x,t;x' ,t') is the Minkowski 
retarded Green's function 



DGret{x;x') 
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— - 



Gret{x;x') = S\x-x') (6) 



and the transverse traceless non-local projection operator is deflned as 



Pij = Sij - ^didj (8) 



where Tij is the conformal coordinate stress tensor and the symbol \pt represents the eval- 
uation at the end of the phase transition. Evaluation of the right hand side is accomplished 
through the power spectrum P{ki,t[,t2) written as 

Tij{t[, h)f*j(t'„ k^)) = (2nf5^'\k, - k,)P{k„ t[, Q [pT%v)\ ' al (9) 

which incorporates spatial translational invariance and a suggestive normalization of bubble 
pressure squared that assumes that the correlator is dominated by disconnected diagrams 
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(i.e. bubble interactions are neglected). Here, pj^^^ is a fiducial energy density of the fluid 
measured by an observer at rest with the fluid at the location of the bubble wall, 7^^. is the 
usual Lorentz contraction factor associated with the fluid velocity Vf behind the bubble wall. 
It is related to the bubble wall velocity by Vf = {v^ — Cg) /(I — VwCg), where Cg is the speed 
of sound and Vy^ is taken to be a constant [42j. We therefore find the gravity wave spectrum 
as 



dp, 



'GW 



1 1 



d\nk (27r)2M2 



rest., 2 
Pf 



dt[dt'2 cos [k{t[ - 4)] [k^P{k, t[,t2)] (10) 



:iii 



Note that the integration is assumed to be over all time, but k^P{k,t[,t2) will have support 
dominantly over a time period At surrounding the time of the PT. As we will check explicitly 
later, this turns out to be reasonable even for relatively longer lived turbulent sources. The 
characteristic sizes governing k^P include the duration of the last part of the PT a^,At < 
(where iJ* is the Hubble expansion rate at the time of the phase transition), the t ime at 
which the PT occurs, and the size of the typical bubble R. Since 

1 1 

R v^a^^At 

if we assume that these are the only important scales in the spectrum and that t^, only enters 
as 

t[-t., (12) 

we can write 

FkAtiiA - Q/At, (4 - Q/At) = eP{k, t[, (13) 

Caveats to this assumption will be discussed in Sec. [5| but the general conclusions there will 
be that this assumption is robust. This gives us 



dp, 



'GW 



dink 



(27r)2M2 
{a.Atf 



(27r)2M2 



rest 
Pf 



1 2 



Vf 



dt[dt2 cos [k{t[ - t'a)] FkAti 



At 



At 



-X14) 



rest 
Pf 



1 2 



Vf 



dq[dq2 cos [kAt{q[ - q'^)] FkAtiq'i, ^2) 



(15) 



where is integrated over (—00, 00) and is dimensionless. Hence, the only conformal sym- 
metry breaking scale in the integrand is At. 

To proceed with the analysis, we need to determine what sets the mass dimension scale 
of At. As is well known [43j, the nucleation rate of the PT bubble per unit volume per unit 
time at temperature T is 



7 = CiT exp 



dt 



/T 



(16) 



where S^^^ corresponds to the appropriate bounce action at finite temperature and Ci is 
assumed to be 0{1). Hence, 



dt 



TrflnT' 

a{t^ 



1 + 



1 dlng,s 
3 dlriT 



dS(^ 
d\nT 



(17) 
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where is the expansion rate at the time t^, of the PT and we have assumed that the entropy 
whose density proportional to g^sT^o^ is conserved where a is the scale factor and g^s counts 



the entropy degrees of freedom. As given by Eq. (A13) of Appendix [A| the completion of 
the PT occurs during a time interval of 



a=i,At oc 



for H^a^^At < 1. Hence, the only Hubble expansion rate dependence in the integrand 

cos [kAt{q[ - gs)] i^fcAt(g'i, ^2) (19) 



is in kAt. 

Explicitly, with the definition 



H7r 



where the U subscript stands for "usual" radiation domination scenario, the integral 



dq'M cos [kAt{q[ - ^2)] ^fcAt(g'i, 52) 



is invariant under the transformation 



At -> At/^ 
k^ki 



(20) 

(21) 

(22) 
(23) 



Hence, with the present assumptions, the spectrum changes under the transformation Hu — )■ 
Hu^ as 



1 1 



(27r)2M2 \a 



p 



rest 2 2 

Pf ^VfVf 



,Atf / rfg;rfg^cos[A;At(g;-g^)]FfcAt(g;,g2) ^ 



(27r)2 M2 V a 
or equivalently 



rest 2 2 

Pf ^VfVf 



a^At/Sf / dq[dq'2CO?,[kAt{q[- q'^)/£]Fk/^t/i{.q'i,q'2) (24) 



(25) 



dpGw{k) _^ 1 dpGw{k/0 



d\nk ^2 d\nk 

This is the main analytic result of this paper. Although this result in some sense has been 
reported in the literature indirectly before (as we survey below), one of the points of this 
study is to examine the robustness of the relationship and to spell out the assumptions 
necessary. 

From this equation, we can thus extract two easy to remember features: 

1. If the peak frequency without quintessence is kp, the new peak frequency is at kp^. This 
is intuitive from recognizing that a smaller bubble length scale results if the expansion 
is faster (corresponding to > 1) since the bubble length scale has a smaller time to 
grow before the PT is completed. 
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2. The amplitude of the spectrum at the peak should decrease as 1/^^ compared to 
case with ^ = 1. The intuition for this result is that less source contributes to the 
gravitational wave at a fixed point because of the past Minkowski null boundary and 
the approximately compact time support of the source. 



Eq. (25) is essentially a classical dimensional analysis coming from the assumptions sum- 



marized by Eqs. (11) and (12), and the result is useful in allowing us to read off how the 
gravity wave spectrum depends on the general assumptions of the Hubble expansion during 
the PT. As we shall see in the literature survey, the details of the spectrum for frequencies 
larger than the peak frequency is difficult to compute and very uncertain. Nonetheless, we 



will show in sections |3j and |5j that our scaling arguments above are robust, and Eq. (25) 
will most likely apply to improved spectra that will be derived by all future more accurate 
computational techniques. 



Note that the amplitude dependence on the fluid energy pj^^^ written in Eq. (24) in a 
very intuitive form is a bit misleading since it naively looks as one can increase the amplitude 
of the gravity wave by increasing this quantity. However, since pj^^^ scales as the radiation 
energy density, when the gravity wave energy density is compared with the radiation energy 
density, one power of it is normalized away. The second power of p^^^^ actually represents 
the clock units with which to measure (a^^At)"^ since when divided by M^, it represents the 
approximate expansion rate squared of a radiation dominated universe. 

Hence, it is not really the absolute magnitude of pj^^^ that is important for increasing the 

measurable gravity wave amplitude but the dimensionless quantity (a*At)^pj^^^/Mp. As we 
will see explicitly in Sec. |4| obtaining a spectrum observable at the LISA experiment will 
require the duration a* At of the PT to be of the same order of magnitude as the expansion 
rate determined by the radiation energy density during the PT. Physically, this effectively 
corresponds to a sufficiently large potential barrier suppressing bubble nucleation such that a 
non-negligible supercooling occurs before the PT completes. Although somewhat tangential 
to the point of our paper, we discuss this issue a bit further in Sec. |4]and Appendix D 



Although our next goal is to utilize Eq. (25) to make predictions for the quintessence 



kination dominated scenario, we will in the next section first check the consistency of our 
result with some of the explicit computations in the literature. There, we will also consider 
turbulence contributions to the gravity wave production and show that even when the tur- 
bulent source is long lived, its dominant contribution will be from the phase transition time 



period, allowing Eq. (25) to remain a good approximation. If the reader is not interested in 



the consistency check, the reader is encouraged to skip to Sec. |4j 



3. SURVEY OF EXPLICIT COMPUTATIONS AND DETAILED ANALYSIS 

In this section, we will survey the literature which computes the gravity wave spectrum 
using both simulations and analytic techniques [3-|5l [TOl [131 [HI [Ml [IH [HI Our aim is 



to show that our scaling assumptions resulting in Eq. (25) are consistent with the existing 
explicit computations. 

References [T3l [HJ [161 [H] find the following generic result regarding GWs generated from 
bubble collision. The spectrum is found to have a rising and a falling shape, where the 
increasing side of the spectrum scales as and the peak position is at a wave vector of 
order of l/i?* {R^ is the typical bubble size at the end of the PT). Physically, one can 
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attribute the scaling law to the compactness of the sources and the spatial homogeneity 
of their distribution. At the same time, one can understand the appearance of the conformal 
symmetry breaking scale as being the only identifiable classical length scale in the 

problem. While there is no well known uncertainty about the rising part of the spectrum 
|18) . there is a large uncertainty associated with the falling part (UV part) of the spectrum. 
Direct simulations find that it scales as fc~^|16]. while the analytic calculations give model- 
dependent results [HI ITS] . 

In addition to bubble collisions, turbulent motion of the fluid during the EWPT will also 
generate GWs [714101 [T9]. One can estimate the turbulence spectrum with either dimensional 
analysis or the velocity correlation function [45j. It is possible that fluid turbulence and 
magnetohydrodynamic (MHD) turbulence can contribute to GW production as much as 
bubble collisions. 

We now consider the details below. 



3.1. Simulation for bubble collision 

GW generation from bubble collisions in zero temperature vacuum has been extensively 
studied by Kosowsky et al. [3j. Each bubble is represented as an 0(3,1) kink solution to 
the scalar field (0) equation of motion. The latent energy released in the phase transition 
is turned into the bubble wall's gradient energy (V0)^ and kinetic energy {dt4>Y. Using 
numerical simulations, it is found that the GW spectrum depends primarily on the large 
scale features of the source. Specifically, a model of the source's stress energy tensor which 
ignores the collision region between the two bubbles still results in a GW spectrum that 
is almost identical to the one from a full simulation. This "envelope approximation" when 
applied to many bubbles collisions |1] {N ~ 20 — 200) results in a GW spectrum characterized 
by a peak frequency and a ratio of GW total energy over the released latent energy as 

oomax = 1.6/3, Ecw/E.ac = OM{Hjpf (26) 

where /3 is a parameter controlling the bubble nucleation rate 7 as in 7 oc exp(/3a*t), and 
/3 is related to the PT duration by /3 ~ {a^:At)~^. Upon going through their arguments, 
one can show that H^: (Hubble expansion rate at the time of the PT) in this equation can 
be seen to be a reparameterization of radiation temperature T^, rather than the expansion 



rate of the universe itself.^ Hence, Eq. (26) is consistent with our result that the fraction of 
the energy of the GW is proportional to the PT duration (a*At)^. A statistical approach 
is also considered in where the gravity wave spectrum is computed as an incoherent 
sum over individual bubbles weighted by the bubble size distribution function. In the many 
bubble case, the high frequency part of the spectrum is enhanced by the multiplicity of small 
bubbles. 

Building on the zero temperature work, the finite temperature situation is studied in 
[5j. Finite temperature modifies the Higgs boson's effective potential and introduces a fluid 
dynamical degree of freedom. If the PT produces large latent heat, the bubble wall velocity 



^ Even though their bubble computations are based on zero temperature vacuum, they add in the expansion 
of the universe by hand to scale the energy density appropriately for the physical values today. This is 
the source of the T* dependence. 
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would be supersonic, leading to a process called detonation. As the detonation front expands, 
the fluid that is swept by the front is being compressed and dragged along (see [T] for details 
of computing the single bubble's velocity proflle and temperature proflle). A full simulation 
of two bubble collision is computationally difficult, due to the chaotic fluid motion and the 
wide range of length scales involved. To circumvent this difficulty, "envelope approximation" 
is again applied, resulting in the following total GW energy fraction and peak frequency [5]: 



n.„.,' « 1.1 X lO-V (-^V (^X (^) ,27) 



(3 J \l + aj \0.24 + vlJ \g, J 
/_«5.2xlO-H.(|)(^)(|^)'" (28) 

where is the detonation front's velocity, is the phase transition temperature, and (7* 
is the relativistic degree of freedom during the EWPT. Here a = Pvac/Prad is the ratio of 
vacuum energy density to the radiation energy density, which characterizes the strength of 
the PT. The variable k ~ p fluid / pvac is an efficiency factor quantifying the fraction of the 
available vacuum energy that goes into the kinetic motion of the fluid. Evidently, the GW 
spectrum's amplitude and peak position have the same scaling dependence on /3 as the zero 
temperature case, which is consistent with our scaling result since /3 quantifles l/(a=i,At) oc H 
while here represents the radiation energy density and not the total expansion rate. (Even 
though here corresponds to Hjj in other sections of the paper, we maintain the original 
literature's notation in this survey section to emphasize the non-transparency of the scaling 
relationship in the literature.) Hence, if exotic fluid component contributions to the stress 
tensor increase the expansion rate during the PT keeping the temperature flxed, then the 
scaling can be read off from the existing literature by scaling with /3. In that sense, the 



change in the peak position and amplitude resulting from our Eq. (25) is not particularly 
new. On the other hand, the robustness of this simple scaling relationship for the entire 
observable spectrum and its application to the kination dominated quintessence scenario 
have not been explored before this paper, to our knowledge. Generalizing the two bubble 
collision simulation to many bubbles (A^ ~ 100) [16], the GW spectrum's high frequency 
part is found to be enhanced from k~^'^ to k~^, again due to the small bubble effect at the 
end of phase transition. 



3.2. Analytic calculation for Bubble Collision 

Reference |14) uses the stochasticity of the source to estimate the gravity wave's spectrum. 
Two assumptions are made about the velocity fleld in the bubble collision: 

1. The velocity fleld's distribution is approximately Gaussian, i.e. the four points corre- 
lator (vvvv) is determined by the two point correlator (ff). 

2. The two point correlator {v{x)v{y)) is nonzero if x and y can belong to the same 
bubble. 
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The following spectrum is thereby obtained: 

dink 

(29) 

where s is the ratio of bubble wall thickness to radius and Z/Zm ~ /c/Zcp is a dimensionless 

wave- number. ^ Since the amplitude contains the factor (^^^ > and the frequency depen- 
dence is through a function of Z = kL^, (where L^, is the typical length scale at the end of the 
stirring phase), the spectrum's parametric scaling is consistent with our scaling rule result 
in Eq. 



3.3. Scaling of the Turbulence Generated Gravity Wave Spectrum 

As the thermal bubbles percolate, the fluid within the bubbles collides and generates 
turbulence. Even though a detailed simulation of turbulence's evolution is difficult, some 
statistical features can be derived from dimensional analysis and intuition. Generally, the 
turbulence contains eddies of different sizes, as the larger ones break down to smaller ones, 
energy is also cascaded down to smaller scales. In a fully developed turbulence, the standard 
intuitive assumptions are that the energy cascade rate is a constant over time and is also 
a constant for different scales |E1 US]. For example, if the cascade rate e were not constant 
on different scales, there cannot be a steady state since energy would be building up at a 
particular length scale. However, as far as eddy sizes are concerned, since a largest scale 
L fixed by the bubble size and a smallest scale fixed by the viscosity exist, the energy 
cascade rate cannot be exactly scale invariant. Assuming that a turbulent eddy of size / 
with velocity Vi breaks down in a few turn over times Ti l/vi (which is true by dimensional 
analysis assuming that viscosity plays a negligible role), the energy cascade rate per unit 
mass for non-relativistic eddies is 

e ~ -y- ~ f?//, or equivalently Vi ~ (e/)^/^. (30) 
l/vi 

It is interesting that the assumptions of e being a constant over different length scales and 
viscosity being unimportant until very short length scales fix the velocity spectrum by di- 
mensional analysis. As / decreases (i.e. for smaller eddies), the dissipation effect becomes 
increasingly important, and the energy dissipation rate per unit mass is given by 

u{yvf^u{^j)" (31) 

where v is the kinematic viscosity. If / is decreased to the point that the energy dissipation 
rate equals the energy cascade rate, the turbulence would cease to exist. This scale is the 
Kolmogorov microscale A: 

e ^ zy (^^ j ' ^ A ~ L(Lt;z,/z/)-3/4 ~ L^Rey^'^ (32) 



^ We use k and w interchangeably in the GW spectrum. 
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where the dimensionless Reynolds number = Lvl/v relate the largest scale L and the 
smallest scale A. In the case of EWPT, the Reynolds number is usually on order of 10^'^, 
causing A to be negligible compared to L, indeed. 

Given these two parameters L and Vl, the energy spectrum of a fully developed turbulence 
can be estimated as 

E{k) ~ e2/^A:-^/^ for A < A;-^ < L (33) 
which is used in [5] to estimate a GW spectrum 

' vvl[ — ] . (34) 



p duj \ P J \oJo 

In the formula, the parameter Vq is the typical fluid velocity at the length scale of the largest 
bubble size L, not to be confused with the bubble boundary velocity f , and ~ r^^^ ~ 
/3f ~^fo. This spectrum is valid up to the smallest scale of turbulence, i.e. ~ A. The 
factors (/3)~^ and u/uq in the spectrum indicate that the turbulence generated gravity wave 



spectrum also is consistent with our scaling rule equation Eq. (25 ). Just as for all the previous 
examples, should be viewed as parameterizing the critical temperature T* and not the 
true expansion rate of the scale factor. 

Reference [19] also considers the initial stirring phase and the final decay phase of the 
turbulence. The velocity correlator function of the turbulence is used to find the stress- 
energy correlation functions. It is found that both the MHD and fluid turbulence last for a 
long time after the PT has ended, in contrast with the bubble collision case. The gravity 
wave spectrum is found to be 

^ go \ f 



dlog k 



\9fiii/ \^^rad*/ 



+ / dy—j-Is{K,,y,y) cos [z - y) 

Jl y+ Jy z + ^ \ Vl 

In this spectrum, the subscript denotes the present time, * denotes the end of PT, and 
fin denotes the end of the turbulence. The subscript s in Cg can be either v or b, which 
stands for the fluid turbulence or magnetic turbulence. Cg is a numerical factor and Qg is 
the corresponding source's energy fraction of the total energy density. The wave-vector k 
is rendered dimensionless as K^, = kL^. The variable y is a time variable normalized by 
the largest eddy turn over time tl: y = where tj„ is the beginning time of the stirring 
phase. The finish time of the turbulence is denoted as yfin{k), the k-dependence indicates 
that different modes end at different times. The integral of z 's upper limit yt^p = min[?/jj„, y+ 
^^r^] serves as a cut-off of correlation between sources at different time. Is{K^:,y,y) is the 
normalized dimensionless equal-time correlator of the source. The index 7 is explained below. 

The first line represents an overall normalization, controlled by the amplitude of the 
source. The second line with time integral dy... represents the contribution to GW from 
turbulence during the stirring up time period, i.e. the PT period. The last line with time 
integral f^^'" dy... represents the contribution to GW during the decay of the turbulence. We 
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shall show the evolution of the turbulence's stress energy tensor only depends on one scale, 
i.e. the PT duration At, which is also the largest eddy turn over time t^. In the stirring up 
part, Tl is clearly the only scale. In the free decay part, there might be a new time scale 
controlling the decay, but we shall see there is none. The decay of the turbulence is modeled 

as 

This power law relation with time is a scale free relation. Therefore, both parts of the 
evolution of the turbulence contain at most one scale. On the other hand the correlation 
between sources at different times is also modeled in a scale free way. It is assumed that the 
source with a certain wave vector k at two different times tiand ^2 are uncorrelated if the 
time lapse is larger than a few oscillation time, i.e. 

{fik,t,)fik,t2))=0, if |ti-t2|> y ,Xe~0(l) (37) 

Thus the evolution and correlation of the source contains no other scale than tl. Finally, 



one can schematically put the above Eq. (35) into the following form: 



Psrc X 



dti / dt2P{kU,ti,t2) 

-OO J 



OO '— J — OO 



where P{kL^:, ^) is a scale free formula. Thus, the GW spectrum in the long lasting 
turbulence case is still consistent with our assumptions leading to our scaling rule equation 
Eq. (|25|. 



4. EXAMPLE: KINATION DOMINATION PHASE OF QUINTESSENCE 



In this section we will apply Eq. (25) to the results of [HI [T6] to compute how the gravity 
wave spectrum will shift due to the assumption of the existence of a quintessence kination 
dominated phase [SSHSS]- The class of models that we are interested in can be described by 
the Lagrangian 

1 

C = -(dqf-V{q) + CM-^R (38) 

where the real scalar field q couples only to gravity (described by the Ricci scalar R) and 
the matter sector Cm (which must contain the Standard Model sector) through the minimal 
metric coupling. As discussed for example in |38) . the quintessence energy density scaling 
with a is not strongly constrained during the early universe if one is willing to tune V{q). If 
the kinetic energy dominates, the phase of the q field is said to be in a kination dominated 
phase, and the energy density behaves as 

P, = + V{q) cc (39) 

with an equation of state w = 1. Starting from this phase, when the kinetic energy has 
decayed away, P/ p can behave as w ~ — 1 equation of state fluid during cosmological periods 
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for which we have empirical evidence for the existence of dark energy. If V{q) participates 
such that it gently pushes q to compensate for the Hubble friction, then pg can decrease less 
quickly than a~^. Instead of focusing on the details of the finely tuned V{q) that can realize 
different scaling behavior with a, we will simply parameterize the quintessence energy decay 



as^ 



Pg oc a"" (40) 

with n G {4, 5, 6}. The results can also be used to understand situations in which the faster 
expansion rate does not arise from quintessence but rather other exotic fluid components 
contributing to the stress-energy tensor [36j. 

A large family of quintessence models can be described by two parameters (n, rj) char- 
acterizing the energy density's scaling behavior (defined above) with Hubble expansion and 
the relative relic density at the time of BBN, respectively. The latter is defined (e.g. ^7 1 138 ] ) 
as 

r, - ^^4r^. (41) 

Current projection of the collider sensitivity to rj if dark matter is an MSSM LSP is 10~^ 
for the LHC (upper bound only) and 10~^ for the ILC ^7\. As we will now see, the gravity 
wave probe sensitivity is even greater. 

We take the gravity wave spectrum from the recent literature and compute the shifted 
spectrum due to quintessence. In Fig. [1} the plots in the left column are based on the formula 
from numerical simulations jlG]. Starting from the nMSSM Higgs model with superpotential 



WnMSSM = XSHi ■ H2 —S + Wmssm (42) 

the result reported in [13j involved a scan over the model parameter space looking for regions 
that give rise to a strong first order PT. One set of FT parameters consists of {a = 0.2, = 
70 GeV, = 30}, taken from set (6) in Table I in (TB]. We shall use this set of parameters 

for all the plots in Fig. [T| To make the plots on the left column of Fig. [T| we use 

n (f)-n (^ + b)f'f'' u-,) 

ncwif) - (43) 



^ See appendix [c| for a formal mapping to potential. 



13 




Figure 1: These plots illustrate the gravitational wave spectrum change due to the effects of a 
kination phase of quintessence. The phase transition parameters are chosen as a = 0.2, T* = 
70 Ge\ , P / Hjj = 30. On the left are plots based on the simulation of Konstandin et al. [16j and 
on the right are plots based on that of Caprini et al. |14| . From the top row downward we have 
n = 4, 5, 6. In each plot, we have the detector sensitivity lines, the original spectrum line without 
quintessence and the lines for rj = 10^^,10~'^,1. 
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where 



(li](£)(^)(i^y" (44) 



/ = 16.5 X lO^^mHz 
^ 0.62 



(45) 

^ 1.67 X 10-.-^.^ r^v r^v r^v' (46) 



+ y/^^ + 2a/3 

^6 « = 7^ 47 

1 + a 

/X 1 r 4 /3a, 

All the variables we do not explicitly discuss below have been defined in an earlier part of 
this paper. The quantities with tildes correspond to the quantities evaluated at the peak of 
the GW spectrum, while those with a star subscript refer to the values defined at the time 
of the phase transition. The parameters a and h in the above first formula correspond to 
the absolute values of the slopes of the increasing and decreasing regions of the gravity wave 
spectrum, and they are fit from the numerical simulations to be a = 3 and 6 = 1 [16]. The 



modified spectra are plotted by shifting the original spectrum according to Eq. (25), where 
^ is computed for each curve from the corresponding (n, rf) as 



Prad{tBBN)V +Prad{tEWPT) -j^'i ^ 

\ ^ ^ , , where pU9. T) = g-T' (49) 

\ PradV'EWPT) JU 



To show that our scaling of the spectra can be applied independently of the details of the 
PT computation, we plot on the right column of Fig. [T] spectra based on analytic estimation 

of m 



^rad^ ?^ ^ 7TT2 / „ x 4.8 



dink 271^ \g J \n*\(3j s 



-d/ ^ l+(£ + £ 



^kin _ 4 (svff 



(50) 



(51) 



^kd 3 1 - {svjr 

s = Cs/vb (52) 

Vf = (vb- Cs)/ (1 - VbCs) (53) 

Qradh'^ = 4.15 X 10"^ (54) 

where Z = kvb/{al3), = 3.8, c,, = a/1/3, g^ = 3.75, g^ = 106.75, and detonation front 
velocity Vf, is related to a as usual. As mentioned above, the PT parameters are chosen to 
be the same as those for the left column plots. 
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The three rows correspond to the parameters n = 4, 5, 6 respectively, and the four curves 
within each plot from top to bottom correspond to rj = 0, 10~^, 10~^, 1 respectively. Ac- 
cording to the bottom row plots, if the underlying Higgs physics is determined sufficiently 
accurately at colliders, a measurement of gravity wave spectrum by the projected BBO ex- 
periment matching the properties of the Higgs model would rule out any kination dominated 
scenario explanation of 0{1%) discrepancy between collider determination of thermal relic 
density and cosmological measurements. More specifically, for 77 -C 1, the discrepancy caused 
by the kination phase can be expressed as 

lO^rj ( ^ ) (55) 



'VlOOGeV 

for a dark matter particle of mass of order 100 GeV |37| . As expected, it is clear from 
the Fig. [1] that GWs have much stronger sensitivity to 7] and n than collider /dark matter 
combination of measurements. Since the kination scenario effectively allows a large (up 
to 10^) boost factor reconciling enhanced galactic annihilations (such as those relevant for 
PAMELA |40| ) with dark matter abundance [38], such measurements from BBO can strongly 
constrain scenarios reconciling collider physics, cosmological DM abundance, and indirect 
dark matter signals. Even more optimistically, if one can obtain the peak position (and the 
right column plot happens to give the correct picture), then one may be able to measure 
both 7] and n, which can be overconstrained by the possible dark matter data which also is 
sensitive to 77 and n. 

Given that the parameters used for Fig. [T] does not lead to a GW spectrum that LISA 
can measure, one might ask "What kind of underlying scalar sector parameters will lead to 

„ ....... . us.r .... ... i^rr _ 

in front of Eq. pA), one simply would have to make this large. However, as we discussed 
at the end of Sec. 2, the duration a^At of the phase transition uses pj^^^ as a clock and 
the two quantities are not independent. For example, in the absence of exotic fluid element 
like quintessence, (a^At)"^ scales inverse proportionally with p^^^^. Hence, the combination 

Pj^^^a*At does not increase even if one increases p^^^^ by itself. Hence, to make the relevant 
combination large, one actually needs a more difficult to achieve model building ingredient 
of making dS^^^ / dT as small as possible. Typically, this requires small thermal corrections to 
the relevant scalar field near the critical temperature. Keeping other parameters fixed, this 
corresponds to a large Tc (this is the temperature at which the PT begins and not when most 
of the PT completes). Qualitatively, this corresponds to a situation in which the system is 
closer to having supercooling rather than not. A semi-quantitative discussion within the 
context of a toy model is presented in Appendix D 

Although we provide no underlying physics model, we plot in Fig. |2] a hypothetical spec- 
trum that would be generated by a "long" PT to demonstrate the scaling effect on a gravity 
wave spectrum measurable by LISA. Because we are setting /3/if = 1 for Fig.|2| this is at the 
edge of the validity of the scaling relationship which assumed that a^H^At < 1. On the other 
hand, the qualitative behavior of the scaling relationship is accurate. Note that a typical 
beyond the SM Higgs sector extended by singlets go through many PTs before the EWPT, 
some of which may lead to the spectrum as shown in Fig. [2j The effect of superposing such 
spectra is beyond the scope of this work. 
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Figure 2: This figure is an analog of Fig. [T] except with the phenomenological parameters of the 
phase transition tuned to give a larger amplitude. The parameters are 13/ Hu = 1, a = 1, T,,, = 70 
GeV. Physically, this corresponds to a phase transition with a long duration {(3/Hu is not large) 
and a large ratio of vacuum energy density to radiation energy density (a = 1 instead of a = 0.2). 



5. CAVEATS 



In this section, we explore the limitations of the assumptions leading to the scaling 
Eq. (25). The most important assumption is that there is only one relevant dimension- 
ful scale, the Hubble expansion rate, during the generation of GWs. Therefore, the existence 
of any other conformal symmetry breaking physics relevant within the kinematic regime of 
our interest may invalidate our argument. The possible scales include the bubble wall thick- 
ness, the correction to the energy density due to bubble interactions, and the dissipation 
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scale associated with turbulence. 

The bubble wall thickness is usually defined near the phase boundary where the energy 
is concentrated. In the vacuum bubble case, the most prevalently nucleated Higgs bubble 
solutions are believed to be described by 0(3, 1) symmetric solutions (f){t,r) which are func- 
tions of — r"^. In the lab frame the bubble wall's thickness will contract towards zero as the 
bubble wall velocity approaches the speed of light. Thus the bubble wall thickness is very 
small compared to the bubble radius. In the thermal bubble case, the energy is modeled 
to be contained in the fiuid's kinetic energy rather than in the Higgs field's profile. Here 
the bubble wall is defined by the portion of fiuid having large enough velocity to give an 
appreciable contribution to GW. Steinhardt [1] showed that there is a large class of solutions 
in which the fiuid's speed only depends on parameter r/t, i.e. as the bubble expands the ve- 
locity profile would expand accordingly. Therefore the bubble wall thickness is proportional 
to the bubble size, and in particular, it is not a scale independent of the Hubble rate. 

Next we would like to consider the effect of small scale physics in the bubble collision 
case. We will show below that the dimensionless gravitational wave spectrum 

eP{kMM) (56) 

in the k range of interest for BBO and LISA depends appreciably only on the dimensionless 
combinations 

where t^, is the time of the PT, is the temperature at the time of the PT which is assumed 
to depend only on the short distance physics parameters M as far as the conformal symmetry 
breaking parameters are concerned, and a* At oc is the duration of the PT proportional 

to the inverse Hubble expansion rate at the time of the phase transition. Since the last of 
these parameters do not scale with the Hubble expansion rate, the conclusion of our paper 



is robust if Eq. (57) are the only relevant conformal symmetry breaking parameters. 

To show this, we start with the effective classical description of the stochastic process of 
bubble creation and stress tensor evolution. Suppose we make the reasonable approximation 
that scalar field evolves classically except with a stochastic source term that can create 
bubbles. Hence, we write the equation of motion for the scalar sector as 

dV d 

^ d(j)* d(j)* ^ ' ' 

where Cj is a short distance physics governed non-derivative interaction Lagrangian, is a 
scalar sector potential, and 

J^{x)= I d'^y7i{y)ji{x - y) (59) 



is a stochastic distribution which accounts for thermal /quantum fiuctuation induced bubble 
creation. The stochastic function Tc{y) is the number of bubbles created per unit time per 
unit volume at a small cell volume^ centered at y and ji{x — y) is a non-stochastic fixed 



^ Bubble nucleation from a classical "appearance" point of view is a collective field process typical of soliton 
creation since a bubble of radius larger than a critical radius is required to appear for it to expand 
classically. 
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function^ representing the effective classical source for a single bubble. Let P(vr) be the 
probability of obtaining a particular function vr for a single realization of the universe. It is 
in principle fixed by the path integral computation of the bubble nucleation, but the details 
will not be important for the scaling arguments that we present if the conditions we later 
discuss are satisfied. The ensemble average {Tij{t[,x)Tij{t2,y)) can be written as 

{Tij{t[,Xi)Tij{t'2,X2)) = J VnP^n) {Tij{t[,Xi)\^Tij{t'2,X2)\n} ■ (60) 



where Tij (a functional of (p) is implicitly dependent on vr through Eqs. (58) and (59) 



(f)i{x) = (f)i[n]{x). (61) 

Note that the functions vr that generate non-vanishing inhomogeneities are the only im- 
portant contributions, even if they may not be at the peaks of the probability functional 
Pin). 

Now, we impose 3 conditions which will limit the number of conformal breaking parame- 
ters that the system depends on: 

1. The value of {Tij{t'i, Xi)Tij{t2, X2)) needs to be known only in the interval t^, < t'^ < 
+ At (where a^At < 1/H^) to obtain an order of magnitude accurate gravity wave 
spectrum since it falls off rapidly outside of that time interval. As shown in Appendix 
\K\ the effective duration of the PT is 

a, At ~ ^ ^ (m/T,(M)) (62) 

where F <C 1 is a dimensionless function of short distance physics conformal symmetry 
breaking parameters M and the critical temperature T^. To make the argument of the 
dimensionless function dimensionless, M has been scaled by the critical temperature 
which itself is assumed to be a function of M only and not H^. The dominance of 
the correlator within this short window of time is a reasonable generic assumption since 
the completion of the PT at + At corresponds to a homogeneous and isotropic scalar 
field phase which cannot source GWs. On the other hand, the plasma inhomogeneities 
can grow (albeit slowly because of the relativistic pressure) through gravitational clus- 
tering. As a first guess, this effect should be most pronounced on short distance physics 
scales (characterized by M) for which the inhomogeneities leading to gravitational po- 
tential would be the largest. If this is true we can neglect these complications as we 
discuss further below. Another caveat is that T* is not necessarily independent of H^, 
since a very large H^, can lead to the decoupling of a particular relativistic species if 
that species has interaction rate T < H^. In that case, the temperature of that species, 
call it Ti, can evolve differently from the rest of the thermal plasma. If a finite Higgs 
VEV can lead to interactions with this species at the time t^, of the PT, one can have 
in addition to T* another temperature Ti{t^) which now does depend on H^,. Note 
that in such situations, T^, and Ti{tt) can be coupled and the dependence on H^, enters 
through Ti(t*). This is highly model dependent and as long as the number of species 



^ Strictly speaking, it is a distribution since it may contain Dirac delta functions. 
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decoupling as a function H^, are not larger than 0(10), its effect will be less than or- 
der 10%. Leaving further investigation of these effects to a future work and staying 
consistent with state of the art computations in the literature, we will not discuss the 
possible breakdown of this assumption further here. 

2. The function 7r(?/) is parametrically dependent on {M, T(t), a^At} where T(t) is the 
time dependent mean temperature of the relativistic fluid. Note that this is certainly 
true for the mean bubble nucleation rate per unit spacetime volume 7(t) . 

3. The functional -P[7r] is parametrically dependent on {M ,T{t),a^,/S.t}. From the com- 
putation of the saddle point approximation of 7(t), this is reasonable. 

These conditions will now be used for a dimensional analysis aimed at checking whether 
there are any physically relevant conformal symmetry breaking scales that we missed. We 
will in effect be reanalyzing some of the steps in Sec. |2] to obtain a sense of what kind of 



details can invalidate the scaling result of Eq. (25). 

First, let's go to the Fourier basis useful for computing the gravity wave spectrum. Because 
of the FRW background assumption, the correlators should have no preferred spatial position 
or direction. Hence, we must have the form 



(Pki 



(T,,(t;,f)T,,(4,y)) = [pTH/f\ < J j^,e^'^-'-^^P{k,AA) (63) 

where the prefactor is guessed from the physical interpretation of the stress energy tensor and 
is defined to be independent of spacetime (these quantities should be expressible in terms of 
M and T^{M)). This guess of the prefactor will be seen to be important below since it gives 
us an argument for the order of magnitude of k^P. Given the assumptions of the previous 
paragraph, we define the dimensionless function F2: 

k^P{k,t[,t'^) = F2{k,T„ At,t[ - U,t'^ - t,, M) (64) 



where we have used Eq. (62), Eq. (17), and assumption number 1 in the previous paragraph. 
Note that the beginning of PT time only appears explicitly in the combination ti—t^,. Given 
that F2 is a dimensionless function, we can make F2 out of the dimensionless combinations 
shown in Table [l} Now, for EWPT occurring at a temperature of order of 10^ GeV, LISA 
and BBO are sensitive to wave vectors at the time of the phase transition in the range 

^|10-'M0-',(M^)-"(^)'"(^)oeV ,5, 

(corresponding to 10^^ to 10^ Hz range detector frequencies). Since we expect F2 ~ 0{1) as 
noted above, for unsuppressed dependence on these parameters, we must have 0(1) para- 
metric possibilities for these dimensionless parameters on which F2 depends. For example, 
for Mia^At ~ 10^^ to produce an 0{1) number consistent with F2, we could have a term 
proportional to 

{Mia,At)-^^^^ (66) 

in which case any scaling associated with At will lead to a suppressed change in F2. Another 
possibility is an extreme fine tuning of parameters 

1/{W-^^ Mia, At). (67) 
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quantity 


range 


quantity 


range 






T*ft' —t*)a* 




kAt 


[10-1,103] 


T, 


0(1) 




~ kAt 


At 

t,;-t. 


0(1) 


k/at. 


[10-15,10-9] 


Mi a* At 


O(10i4) 




O(10i4) 




O(10i4) 



Table I: Dimensionless parameters with a*At ~ 10^^/ff* lOi^GeV^i and T* Mi ~ 100 GeV. 
Since we expect F2 ~ C(l), for unsuppressed dependence on these parameters, we must have 0(1) 
parametric possibihties for these dimensionless parameters on which F2 depends if there is no fine 
tuned dimensionless parameters in the theory. 



This latter possibility cannot be excluded based on dimensional analysis of the form here, 
and our scaling analysis can break down if there are extremely large or small parameters. 
The extreme level of fine tuning must be at least at the level of 1 part in 10i° even when we 
allow for dynamically generated dimensionless numbers of order 10^. Finally, we can divide 
those quantities in the table that are large with each other to give an 0{1) number. Hence, 
from the table, barring an unlikely dynamically generated fine tuning that we discussed, we 
conclude that over the frequency range of interest for gravitational wave detectors 

— * — * 

Since the last factor M /T^{M) is independent of iJ* we arrive at our robust approximation 
that the dimensionless spectrum k^P{k) only depends on the form given in Eq. (13). 



Next, we consider another possible scale, the turbulence's microscale A , which is related 
to the largest scale as A = {Re)~^^^L ~ 10-i*^L. In the short-lasting source's model [5, .lOj , 
this scale serves as a cut-off to the gravity wave spectrum. Since this scale is 10^^ higher 
than the peak scale, we can follow the same reasoning about the short distance scale physics 
to show that it would not be significant for the GW detector LISA or BBO. In the recent 
work [19j where long-lasting source is considered, the microscale can affect the duration of 
the turbulence's free decay. This duration of turbulence may appear as an independent 
time-scale, but as we show in Appendix B, the decay duration is sufficiently long that the 
turbulence has depleted most of its energy towards the end of this duration, and the exact 
ending time, to an excellent approximation, does not enter the GW spectrum. 

Another issue that we did not address are classical scaling violations coming from quan- 
tum radiative corrections. If the C, scaling is many orders of magnitude, the anomalous 
dimension effect may give (depending on what renormalization prescription is chosen) an 
(9(1) correction to the results presented here. However, in that case, GWs from first or- 
der PTs are much less likely to be measurable and therefore are unlikely to be of practical 
interest. 

Before we conclude the caveat section, we would like to comment on the velocity depen- 



dence of the GW spectrum. In [5J, the GW spectrum's magnitude is given by Eq. (27). The 
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GW spectrum depends on the velocity Vu, implicitly through a as for example in 



K 



a 



1 + a 



0.24 + v^, 



(69) 



since Vu, and k are functions of a as in Eq.(47) and (48). In the weak detonation limit of 



a — )■ 0, we have a/1/3 and k — )■ |« , causing the prefactor to scale as a^. In 

Caprini et al. [Tl], the velocity dependence goes as 



(svf) 



2\4: 



(70) 



which in the weak detonation limit scales as (since Vf ^ A/3a;/2, s — ?■ 1 — y/2a). For 
other region of parameters, as is shown in Fig. (13b) of [14j, the two approaches can give peak 
amplitudes that differ by one order of magnitude. It is clear that the velocity dependence 
of the GW spectrum is both uncertain and can be of high polynomial power, making the 
current numerical uncertainty in the bubble wall velocity (see e.g. ^5]) a significant source 
of overall GW spectrum uncertainty relevant for assessing the measurability of the GW. It is 
comforting however to know that our scaling rule is largely independent of this uncertainty. 



6. SUMMARY 



In this paper, we have presented an analytic transformation rule Eq. (25) that is useful 
for understanding how the gravitational wave spectrum generated through an electroweak 
scale first order phase transition at a fixed temperature T^, would change if the expansion 
rate of the universe during the phase transition were different from that inferred from the 
assumption of pure radiation domination. We have explored the remarkable robustness of 
the scaling relationship with respect to many computational uncertainties in the gravity 
wave spectrum. 

We apply this transformation rule to the example of a universe having a quintessential 
kination dominated phase and find as expected a strong sensitivity to the single phenomeno- 
logical parameter controlling this scenario. In principle, this scaling relationship together 
with dark matter properties measured by colliders can be used to overconstrain this single 
phenomenological parameter. Unfortunately, if the current technology of gravity wave com- 
putations is correct, then we find that any measurement of the gravity wave spectrum at the 
level of the projected BBO sensitivity can rule out any appreciable boost factors relevant 
for reconciliations between various sets of data such as that between colliders and cosmology 
and/or indirect detection (such as that relevant [38j for PAMELA data [40j). 

Nonetheless, using the results of this work, any future gravity wave detection experiments 
measuring phase transition induced gravity waves can understand their measurement's sen- 
sitivity to the expansion rate of the universe. It would indeed be exciting to have an ob- 
servational anchor on the expansion rate of the universe when the universe is as hot as 100 
GeV, just as isotope abundance measurements allow us to have an observational anchor on 
the expansion rate at a temperature of 1 MeV in the context of big bang nucleosynthesis. 
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Appendix A: Bubbles Filling Space 



Here we review a well known argument IIB] about how first order PT bubbles fill space. 
For this section, we will use the metric parameterization ds^ = dt^ — a^{t)\dx\^. Let '-jit) 
denote the probability per volume per time of bubble fomation. Assuming that the bubble 
wall is not accelerating in the locally inertial frame (which is what is typically done in the 
literature when the bubble wall velocity is taken to be a constant), we have 

dx I , . , 

— oc -. Al 
dr a 

Hence, we have 

1 dr K 
^ = - (A2) 

where the constant K is to be determined by a boundary condition. Setting the boundary 
condition 

\ /AON 
flj-^li, = Vyj, (A3) 

we find 



K = -^^^ (A4) 
dr Vw 



and 

(A5) 



dt a 

Hence, a bubble nucleated at time ti fills a comoving volume 



47r 

Vz{U,t) = — 



dt' 



-\ 3 



(A6) 



Next, we compute P, the probability that a point in comoving space is in false vacuum. 
The probability that at time t -\- dt the vacuum at a point is still in false vacuum given that 
it is in false vacuum at time t is 

P{t + dt) = P{t)[l - P^] (AT) 

where P^ is the probability of nucleating a bubble within the past causal cone surface volume 
of thickness dt with the causal signal propagation speed given by v^, since it is the bubble 
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wall that needs to reach the given point in consideration. Since the relevant surface volume 
of the causal cone can be easily computed to be 



aV3(t„t)cit, 

we can multiply this by 7 to find Pc to arrive at 

dP 

Solving for P, we find 



P{t) = P{U)exp {- j^^ 



dt'^{t')V^{U,t')a\r 



(A8) 



(A9) 



(AlO) 



Using Eq. (16), we find 
P{t) = P{ti) exp ^-Ci J^^ dt' exp ^ + 

When t < tc where 

~ ^* + 



{t'-QH, rfg(3) 

^3 dlnT 



din ^(3) 

I Ri dings dlnT 
' 3 dlnT 



-\ -1 



\t* 



T\t')Vs{U,t')a%t') ] . 
(All) 

(A12) 



Si'^^ /T is large such that the exponential suppression in the integrand makes the integral in 
front of Ci negligible. After that time scale, probability of not being in the false vacuum is 
0{1), and the PT is assumed to be completed. Hence, we can conclude that the duration of 
the PT scales as 

1 



At 



proper 



In terms of conformal time, when i^^Atproper ^ 1, we have 



^^conformal ~ ^^proper/a(t*) oc 



1 



(A13) 



(AM) 



The variable ^^^^conformal corresponds to the variable At in Eq. (18). 



Appendix B: The decay duration of turbulence 

In the GW spectrum formula Eq. (35) in reference [19j , the free decay part of turbulence 
gives a contribution that proportional to the integral 



Here, y and z are dimensionless time variables: e.g. y = (t—tin)/TL where tj„ is the beginning 
time of the stirring phase and tl is the largest eddy turn over time. We now show that the 
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ending time of turbulence as a function of the scale length yfin{k) is very large for the peak 
position {yfin{kpeak) ~ 10^) and therefore the exact ending time is irrelevant to GW spectrum 
at the peak position. We shall consider the behavior of the integrand as y ^ yjin , and take 
7=K = 2/7 for concreteness: 

dy^IsiK.^y^y)]^ ^cos^^(.-,)j~y^ dyy-^-OiyjtJ (B2) 

where we have used the following estimation 

^ ^ ^ 1^-105 turbulence 

^ *' ' \y^^ fluid turbulence 

ytop -y = min[?//i„ - y, -^] ~ vl/K^ ~ 0(1) 

VTA* 

Therefore the GW's relative dependence on yfi^ is as weak as yj^^ . For the ending time of 
turbulence at the stirring scale K^: = kL^, ~ 1 , one can use Eq. (72) of |19] 

14 ^ 

yf^n - ^ - { , for > 0.07 (B3) 



tl \ K* 

which shows that yfin{K^: = 1) ~ 10^. If one take yfin — )■ oo in the above integral, the error 
introduced is on the order of 0(10^^^), very small compared to the integral itself which is 
at least 0(1). 



Appendix C: A Formal Map between V{q) and n 

Here we give a formal map between the quintessence potential V{q) and the quintessence 
energy density dilution behavior a~". 

The equation of motion of the quintessence can be written using the scale factor a as a 
time variable assuming that the only other component during the era of interest is radiation 
which dilutes as a~^: 

1 V2VV-(g)/p^o + (^)4 d f V2^Viq)/pm + ifrd{q/Mp) ^ ^ MpV'{q) _^ 

fc 1 ( dq \2 din a 1 fc l ( dq \2 din a / ppn 

The solution to this equation defines a functional 

qs[a,V{q)]. (C2) 

Note that the equation has been normalized to be dimensionless: V is measured in units of 
initial radiation energy density Prq and q is measured in units of Mp. Next, one can solve 
q = qs[a,V{q)] for a[q,V{q)]. This can be put into the energy density scaling equation for 
the quintessence energy density, yielding the equation 

'V2jV{q)/pm + 



PqO / '^0 Nn ^ 1 

PRO^a[q,V{q)Y 2 



^a[q,V{q)]' 



c _ ( dqs [ay{q)]/Mp \2 
^ V d\na ' 



, dq,[a,V{q)]/Mp ^^ V{q) 
din a ^'^='^['''^(^)] + 



(C3) 

which can in principle be solved for V{q). 
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Appendix D: An Estimate of Temperature Dependence of S^^^ 

In this section, we will give a semi-quantitative argument of which parameter choices 
of effective potential governing the PT will lead to an enhanced gravity wave amplitude. 
Consider the high temperature expansion of the effective potential for a single real field in 
the form 

V{<P) + cT')4>' - EcP' + (Dl) 

where //^ < and c is a thermal correction dependent parameter. Equation of motion yields 

= + il/ + cT^)(P - 3E(P^ + X(j)\ (D2) 

where we have estimated V </> ~ Yji^'- solving for we need a characteristic value for (f) 
which we will call 0c- We can set this characteristic value to be between the local maximum 
0„ and the minimum 0* (not the one at the origin) of the effective potential. (Recall that 
the cubic term —Ecj)^ is responsible for there being a bump in the potential giving rise to a 
local maximum.) Explicitly solving 1^' (0) = 0, we find 




and 



Hence, if we make a somewhat arbitrary but reasonable definition for the characteristic value 
to be 

fe^^.(r) + ^-'^';^°'^) . (D5) 

we find 

op 

= 2A (°«> 
independently of the temperature except through E and A which we assume to be dominated 
by the non-thermal contribution [E coefficients that rely on thermal corrections do not yield 
strong phase transitions typically anyway) . This yields the length scale associated with the 
bubble action to be 

L = , ^ ^ (D7) 

where we have displayed our assumption of //^ < manifestly. Hence, we have 

S'^'^ ^ i^^((|)^ + V^(0c)) (D8) 

(D9) 
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Computing the temperature derivative at the critical temperature of 



T = Tc = ^ Z- ' (Dll^ 



(which is obtained by setting V{(f)^,) = and solving for T) we find 



T=n ~ 50^V2 + A|/i2|/E2. (D12) 



dT A 
Hence, we conclude 

F ~ -|Li:jfV2 + A|,^|/E^. (D13) 

To make this order unity (appropriate for an enhanced gravity wave amplitude), we can 
consider y/c <^ 1. Assuming A ~ C(l), we find 

c < 10-^ (DM) 

This leads to a PT temperature of 

T.>Tevf^L)vM!±A 

VlOOGeVy ^ ^ 

corresponding to a high temperature PT of a weakly coupled scalars. 
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